Separability of very noisy mixed states and implications for NMR quantum computing 
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We give a constructive proof that all mixed states of A'' 
qubits in a sufficiently small neighborhood of the maximally 
mixed state are separable. The construction provides an ex- 
plicit representation of any such state as a mixture of prod- 
uct states. We give upper and lower bounds on the size of 
the neighborhood, which show that its extent decreases ex- 
ponentially with the number of qubits. We also discuss the 
implications of the bounds for NMR quantum computing. 

PACS numbers: 03.67.-a, 03.67.Lx, 76.60.-k, 89.80. +h 



In this Letter we investigate the structure of the space 
of density matrices of N spin- 1/2 particles (qubits). In 
particular, we consider density matrices that are close to 
the maximally mixed density matrix and ask whether or 
not they are separable, i.e., whether they can be writ- 
ten as mixtures of direct products. One might imagine 
that the issue is straightforward, in that the maximally 
mixed state seems to be very far from the boundary be- 
tween separable and nonseparable states. It might be the 
case, however, that the maximally mixed density matrix 
is surrounded by separable matrices, but that these sep- 
arable density matrices lie in a low-dimensional subspace 
within the space of all density matrices. By leaving this 
subspace, even infinitesimally, one could reach entangled 
density matrices. 

In Q this problem is addressed by an existence proof; 
namely, it is shown that there exists a sufficiently small 
neighborhood of the maximally mixed density matrix in- 
side which all density matrices are separable. In a 
lower bound on the size of the neighborhood is given. 
Here we go further by giving a constructive proof that 
provides an explicit representation of any state suffi- 
ciently close to the maximally mixed one as a mixture of 
product states. We give an upper and a much improved 
lower bound on the size of the neighborhood, which show 
that it decreases exponentially with the number of qubits. 

Our results have immediate implications for present re- 
search that makes use of high-temperature nuclear mag- 
netic resonance (NMR) for quantum information process- 
ing and quantum computation [|3|-p^. Since the first 
proposals to use NMR for quantum computation, there 
has been surprise about the apparent ability to perform 



quantum computations in room-temperature thermal en- 
sembles. It has been a puzzle how these thermal states, 
which are very close to the maximally mixed state, could 
correspond to truly entangled states ||l^. The bounds 
we calculate show that all states so far used in NMR 
to simulate quantum computations or other quantum- 
information protocols are separable. This is not to say 
that NMR techniques are incapable of producing entan- 
gled states, in principle. Increasing the number of corre- 
lated spins might eventually lead to nonseparable states, 
but this question is left open by the bounds derived in 
this paper. 

We consider arbitrary density matrices for N qubits, 
written as 



= (1 - e)Md + epi 



(1) 



where d = 2^ is the Hilbert-space dimension for N 
qubits, Md ~ Id/d is the maximally mixed density ma- 
trix (Irf is the identity matrix in d dimensions), and p\ is 
an arbitrary density matrix. Any density matrix can be 
written in the form (|l|). We show that for e sufficiently 
small, all density matrices of the form (|l|) are separable. 
We define two kinds of representations of p^ in terms of 
product states, which provide candidates for ensemble 
decompositions of pe as a mixture of product states. By 
considering these candidate decompositions, we derive an 
explicit lower bound on the size of the neighborhood of 
separable states. We conclude by establishing an explicit 
upper bound on the size of the neighborhood. 

Our approach is to represent an arbitrary density ma- 
trix in an ovcrcomplete matrix basis, each basis element 
of which is a pure direct-product density matrix. If all 
the coefficients of a density matrix in this representation 
are nonnegative, the coefficients can be considered to rep- 
resent probabilities, and the density matrix is separable, 
as it is then a mixture of direct products. 

All of our representations arise ultimately from ex- 
panding a density matrix for IN qubits in terms of direct 
products of Pauli matrices: 

1 



(2) 



Here and throughout we sum over repeated indices: 
Greek indices run over the values 0,1,2,3, and Latin in- 
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dices take on the values 1,2,3. The matrix ctq = I2 is 
the two-dimensional identity matrix, and the matrices 
(Ti, i — 1,2,3, are the Pauli matrices. The (real) expan- 
sion coefficients in Eq. are given by 



= tr(p cr„^ 



(3) 



Normalization requires that C0...0 — 1- Since the eigen- 
values of the Pauli matrices are ±1, the expansion coef- 
ficients satisfy 



1 — Cai...ajY — 1 



(4) 



To be concrete, we consider first the case of two qubits. 
For each qubit we introduce six pure density matrices, 
Pi = ^(l2 + CTi) and Pi = i(l2 — <Ji). A convenient dis- 
crete overcomplete basis for discussing separability con- 
sists of the 36 direct-product projectors, each of which is 
a pure direct-product density matrix: Pi (g) Pj, Pi (g) Pj, 
Pi (g) Pj , Pi (g) Pj . Any density matrix of two qubits can 
be expanded in this basis, but since the basis is over- 
complete, the representation is not unique. We make a 
specific choice, as follows. Noting that (T; ^ Pi — Pi and 
I2 = Pi + Pi, we can write I2 — coiiPi + Pi), where 
uji — 1/3, i — 1,2,3. With these results we can convert 
the Pauli representation into the form 

P = ^ [(wjWj + CioUJj + LUiCoj + Cij)Pi (g) Pj 



+ {uJiUJj - CioUJj + LUiCoj ~ Cij)Pi (g) Pj 

+ {uiujj + CiQUj - uJiCoj - Cij)Pi (g) Pj 

+ {uJiLUj - C^qUJj - UJiCOj + Cij)Pi (g) Pj] 



(5) 



If the coefficient of each of the 36 basis elements is non- 
negative, the density matrix is separable. We note that 
when the maximally mixed density matrix for two qubits, 
M4 = ■|l2<g'l2, is represented as in Eq. (||), the coefficient 
of each of the basis matrices is 1/36. 

Consider now an arbitrary entangled (nonseparable) 
density matrix pi. Since pi is entangled, at least one of 
the coefficients in the representation of pi in the form (j^) 
is negative. Suppose now that pi is mixed with the 
maximally mixed density matrix M4 as in Eq. (Q), i.e., 
Pe — {1 — e)Mi + epi. Although some of the coefficients 
of pi are negative, all of the coefficients of M4 are strictly 
positive. Hence, for e small enough, all the coefficients of 
Pe are nonnegative, making separable. Thus all den- 
sity matrices in a sufficiently small neighborhood of the 
maximally mixed density matrix are separable. 

Furthermore, we can find an explicit bound on e such 
that pe is separable for any pi. To find a bound, we use 
Eq. (^) to bound the coefficients of the basis matrices in a 
representation of pi of the form . The minimum value 
of any of the coefficients is (1/4) (1/9 - 1/3 - 1/3 - 1) = 
— 14/36. Thus all the coefficients of the density matrix 

in the discrete overcomplete basis are nonnegative if 



(1 - e)/36 - 14e/36 > 0, i.e., if e < 1/15. For e < 1/15, 
the representation (^ is an explicit decomposition of pe 
as a mixture of direct products. 

A similar analysis can be carried out for any number of 
qubits. Starting from the Pauli representation (^), we in- 
troduce a discrete product basis, like that for two qubits, 
and define a representation analogous to that in Eq. (||). 
Using Eq. (^ to limit the size of the coefficients in this 
representation, we find an asymptotic lower bound on 
the size of the neighborhood of separable density matri- 
ces that is of order e ~ 1/4^ for N qubits. 

One particularly interesting example is the GHZ state 
[|l2|,0, a state for three qubits whose density matrix is 

Pghz 

= i(|lll) + |222))((lll| + (222|) 

= ^ (I2 <Xi I2 <Xi I2 + I2 (Xi 0-3 (Xi ^3 -I- 0-3 (g) I2 <gi ca 

+ 0-3 (Xi era (g) I2 + (Ji (gi (Ti (gi CTi - CTi (gi 0-2 (gi 0-2 

- (72 (Xl (Tl <gl Cr2 - 0-2 0-2 <8l 0-1 ) ■ (6) 



We now express the maximally mixed density matrix, 
Mg, and the GHZ density matrix in terms of an over- 
complete set of 6"^ — 216 basis matrices analogous to 
the 2-qubit matrices introduced above. We find that 
Ms has coefficient 1/216 for all the basis elements and 
that the smallest coefficient for pcHZ in this basis is 
-(l/8)(26/27) = -26/216. Thus for e < 1/27, the state 



Pe ^ (l ~ e-)Ms + epGHZ 



(7) 



is separable. We return to the GHZ example below. 

We have also considered another overcomplete basis 
for the space of density matrices, a basis labeled by con- 
tinuous parameters. An arbitrary density matrix for N 
qubits can be represented as 



p = dfli - ■ ■ dilN w{ni,. . . , nN)Pni 



Prn : (8) 



where the integral runs over N Bloch spheres and where 
Pfi = ^(l2 + n ■ (j) is the projector onto the pure state 
located at unit vector n on the Bloch sphere. The rep- 
resentation (^) is by no means unique. In a spherical- 
harmonic expansion of w{ni, . . . ,nN), the density ma- 
trix determines only the / = and I = 1 parts; the 
higher-order spherical-harmonic content corresponds to 
the freedom in representing p as a sum of one-dimensional 
product projectors. A separable density matrix is one for 
which there is an expansion such that w{ni, . . . ,n]\f) is 
everywhere nonnegative. 

We can generate a candidate for a separable ensemble 
decomposition of p by considering the unique represen- 
tation of the form (j^) such that w{ni, . . . , n^r) has only 
I = and I = 1 components. We can obtain this unique 
representation by noting that 
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1 3 



(9) 



where ng = 1/3. Inserting this result into the Pauli- 
matrix expansion (||) and using Eq. (^) gives 



w{ni, 



...,nN) 
3 \^ 



(4^) 



X-^tr(^p{l2 + 3ni • (?) (g) • • • (g) (I2 + Sun ■ ff)^ 



(10) 



The maximally mixed density matrix, M2N , has w — 

Let us concentrate on the operator product in the last 
form of Eq. (^0|) . Each operator in the product has eigen- 
values 4 and —2. Thus the most negative eigenvalue of 
the operator product is 4^^^(— 2) = —2^^^^, which im- 
plies that 



w{ni, . . . ,nAr) > 



-)2N-1 



(4^) 



N 



(11) 



Consider now the density matrix Its candidate 

ensemble probability satisfies 



We{ni, . . . ,n7v) 



1 - e 



-I- eiyi(ni, . . .^un) 



> 



(47r)^ 

1-6(1 + 22^-1 



(4^) 



N 



Therefore pe is separable if 
1 

e < 



1 + 22^-1 N^oo 4^ 



(12) 



(13) 



We see again that all density matrices in the neighbor- 
hood of the maximally mixed density matrix are sepa- 
rable, and we obtain a lower bound on the size of the 
separable neighborhood, which for large N is much bet- 
ter than the bound, e < (1 -t- 2^-i)-(^-i), given in g]. 

It is instructive to return to the GHZ state (||) and to 
note that Eq. (HOh gives 



wghz(^i,^2,^3) 



1 



> 



26 



(47r)3 



■ [1 -1- 9(ciC2 + C2C3 -I- C1C3) 

+ 27S1S2S3 COs((pi +(p2+ m)] 

(14) 



Here cj = cos 9j and Sj = sin 9j , and the minimum value 
occurs at ^1 — O2 = ^3 = tt/2 and ipi + ip2 + (ps = tt. 
Equation (|l^) shows that the mixed state (^ is separable 
if e < 1/27, the same bound obtained above. 

This bound is not optimal. To find a bound for a par- 
ticular state, such as the GHZ state, one should expand 



the state in terms of a tailor-made set of direct products, 
instead of a general-purpose set. The continuous set in 
Eq. (||) provides a starting point for developing a more 
efficient representation (treated in an upcoming publica- 
tion) that uses a linearly independent set of 4^ direct 
products of the form Pft^ (g • • • (g) P^j^ . The explicit form 
of this representation is 



P = 



E 



w{ni, . . .,nN)Pfi^ 



Pn. , (15) 



where w{ni, . . . , njv) is given by Eq. (p^, and the sum 
runs over Bloch vectors that lie at the vertices of tetra- 
hedra. Using a representation of this sort matched to the 
GHZ state, one can improve the bound for separability 
of the state (0) to e < 1/(3 + 6^2) ~ 1/11.5. 

In this Letter we have been using e to characterize how 
close a density matrix of the form (^ is to the maximally 
mixed density matrix, Md- An alternative distance mea- 
sure, defined by (5 = ^tr((p — Afd)2), leads to similar 
overall conclusions, for using the representation of p in 
Eq. (p^), one can show that all states with 5 < 1/(2^5)^ 
are separable. 

Up to this point we have been thinking of the num- 
ber of qubits as being fixed, and we have investigated 
the boundary between separability and nonseparability 
as the amount of noise, specified by e, changes. We now 
shift gears, thinking of the qubits as particles with spin 
and asking what happens as the number of particles or 
their dimension changes, while e is held fixed. In gen- 
eral, as we go to more particles or higher spins, we find 
that we can tolerate more mixing with the maximally 
mixed state and still have states that are not separable. 
In other words, for a given e, we can always find states 
of sufficiently large numbers of particles or sufficiently 
high spin for which is nonscparable. We translate this 
result into an upper bound on the size of the separable 
neighborhood around the maximally mixed state. 

Consider now two spin-((i — l)/2 particles, each liv- 
ing in a c?-dimensional Hilbert space. What we have in 
mind is that each of these particles is an aggregate of 
N/2 spin- 1/2 particles (qubits), in which case d = 2^1"^. 
We consider a specific joint density matrix of the two 
particles. 



= (l-e)Af<,2+e|V;)(V| 



(16) 



where \i/j) is a maximally entangled state of the two par- 
ticles, 

|V;) = -^(|l)|l) + |2)|2) + ... + |d)|d)) . (17) 

Now project each particle onto the subspace spanned 
by |1) and |2). The state after projection is 

1 /I - e 

P 



>|1> + |2>|2>)((1|(1| + (2|(2| 
= (l-e')M4 + e'|<A)(<^| , (18) 



A\ .2-1^ + ^(11^ 
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where A = {A/(P)[1 + e{d/2 — 1)] is the normahzation 
factor, 



|0) = i=(|l)|l) + |2)|2)) 

is a maximahy entangled state of two qubits, and 

ed/2 



^ _ 2e/d 



A l + e(d/2- 1) ■ 



(19) 



(20) 



The projected state p is a Werner state [|l7|], a mix- 
ture of the maximally mixed state for two qubits, Af4, 
and the maximally entangled state \<j)). The proportion 
e' of maximally entangled state increases linearly with 
d. Thus, as d increases for fixed e, there is a critical di- 
mension beyond which p becomes entangled. Indeed, the 
Werner state is nonseparable for e' > 1/3 p7| , p^ , which 
is equivalent to d > e^^ — 1. Moreover, since the local 
projections on the two particles cannot create entangle- 
ment from a separable state, we can conclude that the 
state ( p^ ) of N qubits is nonseparable under the same 
conditions, i.e.. if 



e > 



1 



1 



1 + 2^/2 



(21) 



This result establishes an upper bound, scaling like 
2^^^^, on the size of the separable neighborhood around 
the maximally mixed state. 

Our results have implications for attempts to use 
high-temperature NMR techniques to perform quantum 
computations or other quantum-information-processing 
tasks. They imply that NMR experiments performed to 
date have not produced genuinely entangled density ma- 
trices. This is because in current experiments, the pa- 
rameter e, which measures the deviation from the maxi- 
mally mixed state, has a value 10^'°^, much smaller than 
the lower bounds we have found for the radius of the sepa- 
rable neighborhood of the maximally mixed state, for the 
cases of two or three spins used in these experiments. 

Present high-temperature NMR techniques, based on 
synthesizing a pseudopure state in the deviation density 
matrix, imply that e scales like N/2^ as the number of 
qubits increases at constant temperature ||^,^ . With this 
scaling, the state leaves the region where our lower 
bound implies that all states are separable at about 14 
qubits, but it never enters the region where our upper 
bound guarantees that there are entangled states. Thus, 
it is unclear whether present NMR techniques can pro- 
duce entangled states. Different techniques might lead to 
a more favorable scaling behavior for e JTot . 

The results in this Letter suggest that current NMR 
experiments should be considered as simulations of quan- 
tum computations rather than true quantum compu- 
tations, since no entanglement appears in the physical 
states at any stage of the process |20j. We stress, how- 
ever, that we have not given a proof of this conclusion, 



since we would need to analyze further the power of gen- 
eral unitary operations in their action on separable states. 
Much more needs to be understood about what it means 
for a computation to be a "quantum" computation. 
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